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ABSTRACT 


This  is  the  second  part  of  an  investigation  dealing  with  the  global  load- 
displacement  behavior  of  a  center  cracked  panel.  Energy  dissipating  mechanisms 
of  microcrack  formation  and  macrocrack  propagation  are  modeled  separately  and 
in  combination  in  accordance  with  a  predetermined  criterion  providing  the  limits 
for  each  mechanism.  Nonlinear  global  behavior  of  the  cracked  panel  is  exhibited 
on  the  load-displacement  curve  developed  by  incremental  loading  and  macrccrack 
growth  with  or  without  microcrack  damage.  .  ~ 

The  regions  of  microcrack  initiation  to  the  sides  of  the  macrocrack  are  as¬ 
sumed  to  coincide  with  the  locations  of  maximum  strain  energy  density  (AW/AV)max 
condition  similar  to  that  of  Haigh  and  Beltrami  for  yielding.  The  material  ele¬ 
ments  within  these  regions  modeled  by  the  finite  element  grids  undergo  a  decrease 
in  the  local  stiffness  which  reflects  the  severity  of  damage.  The  damage  thres¬ 
hold  for  each  element  corresponds  to  the  elastic  strain  energy  density  at  the 
yield  stress  level  of  a  uniaxial  tensile  test  specimen.  As  the  material  elements 
are  damaged,  the  initial  linearly  elastic  behavior  changes  to  pseudo-linear  elas¬ 
tic  behavior,  under  the  assumption  that  any  energy  dissipated  due  to  microcrack 
generation  is  no  longer  recoverable  upon  unloading.  No  other  dissipative  mecha¬ 
nism  of  microdamage  is  assumed.  The  proposed  Multiple  Damage  model  describes 
the  local  effect  of  microcracking  by  a  finite  element  formulation  that  utilizes 
twenty-four  (24)  discrete  values  of  the  secant  elastic  modulus  of  a  specified  uni¬ 
axial  tensile  true  stress-true  strain  curve.  The  most  severe  level  of  material 
damage  corresponds  to  a  secant  modulus  evaluated  from  the  point  of  ultimate  stress 
and  strain.  This  represents  the  effective  properties  of  the  microcracked  portion 
of  the  medium. 
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The  direction  of  macrocrack  growth  is  also  assumed  to  be  governed  by  the 
Strain  Energy  Density  Theory.  In  contrast  to  microdamage,  the  locations  of 
( AW/ aV )  .j n  determines  macrocrack  advance.  The  increment  of  crack  growth  at  a 
given  load  depends  on  the  critical  magnitude  of  AW/AV  characteristic  of  the  ma¬ 
terial  ahead  of  the  crack  front  and  is  designated  by  (aW/aV)  in  the  absence  of 
damage  at  the  microscopic  level.  Otherwise,  the  local  stiffness  of  the  material 
along  the  path  of  expected  crack  growth  may  also  be  affected  and  should  be  dis- 
tinguished  accordingly  with  a  different  notation  (&W/aV)c.  The  value  of  (AW/AV)c 
is  a  relative  value  which  represents  the  decreased  resistance  to  fracture  indica¬ 
tive  of  the  degree  of  material  damage.  This  is  consistent  with  those  situations 
where  microdamage  ahead  of  the  macrocrack  precedes  the  macrocrack  propagation. 

The  material  is  weakened  in  that  less  energy  is  required  to  advance  the  crack 
front  which  corresponds  to  a  reduction  in  (aW/aV)  . 

The  applied  load  is  considered  as  the  independent  variable  which  increased 
incrementally  while  material  damage  at  both  the  microscopic  and  macroscopic  level 
is  monitored  by  the  model  described  above.  The  stress  and  failure  analysis  is 
repeated  for  each  increment  of  macrocrack  growth.  This  could  be  done  up  to  the 
point  of  fracture  termination  that  is  the  final  total  separation  of  the  panel, 
assuming  that  the  material  parameter  governing  crack  instability  is  known.  The 
results  for  two  different  fracture  toughnesses  are  displayed  graphically  and  dis¬ 
cussed  in  terms  of  the  degree  of  load-displacement  nonlinearity  caused  by  nonuni¬ 
form  crack  growth  rate  and  change  of  local  stiffness  of  the  material  owing  to 
microcracking.  It  is  found  that  the  individual  contribution  of  these  two  effects 
depends  on  the  fracture  toughness  of  the  material  while  the  loading  rate  was  kept 
constant. 
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INTRODUCTION 


The  interaction  of  material  damage  at  the  microscopic  level  accompanied  by 
macrocracking  in  engineering  materials  is  of  fundamental  interest.  A  basic 
understanding  of  this  can  lead  to  an  extension  of  the  admissible  range  of  allow¬ 
able  load  for  structural  components.  This  requires  the  identification  of  the 
successive  failure  modes  with  the  load  history  of  a  given  structural  component 
which  behaves  in  a  ductile  manner.  The  optimization  of  strength  and  ductility 
cannot  be  achieved  without  a  consistent  and  reliable  procedure  that  can  qualita¬ 
tively  and  quantitatively  account  for  the  path  dependent  material  damage  process. 
This  involves  a  description  of  slow  and  fast  crack  growth.  It  is  well-known  that 
the  local  stiffness  of  a  material  is  reduced  when  yielding  takes  place.  Micro¬ 
photography  of  the  yielded  materials  [1]  shows  that  mechanical  imperfections  in 
the  form  of  minute  cracks  or  cavities  are  formed  which,  in  this  investigation, 
will  be  referred  to  as  microdamage.  This  type  of  damage  that  has  been  commonly 
observed  in  regions  close  to  the  macrocrack  front  can  significantly  influence 
the  growth  characteristics  of  the  macrocrack  and  the  load  at  fracture  instability. 

The  classical  continuum  theory  of  plasticity  has  not  had  much  success  to  ex¬ 
plain  the  foregoing  physical  damage  process  for  a  number  of  reasons.  One  of  them 
is  that  the  yield  condition,  say  the  Huber  or  Henky-von  Mises  criterion,  is  cho¬ 
sen  independent  of  the  fracture  initiation  and/or  propagation  condition  of  the 
macrocrack  [2,3].  This  can  introduce  a  large  degree  of  inconsistency  and  arbi¬ 
trariness  into  the  stress  and  failure  analysis.  An  attempt  is  made  here  to  model 
microscopic  material  damage  near  a  macrocrack  by  a  finite  element  procedure  in 
which  the  local  stiffness  at  the  different  locations  change  in  accordance  with 
the  strain  energy  density  criterion  that  can  uniquely  predict  microdamage  or 
yielding  as  well  as  the  growth  of  the  macrocrack.  It  has  been  established  pre- 
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viously  [4,5]  that  the  relative  maxi mums  of  AW/aV  correspond  to  yielding  and 
the  relative  minimums  to  fracture.  In  this  way,  a  consistent  procedure  can  be 
developed  to  model  material  damage  by  a  combination  of  fracture  and  yielding, 
a  situation  often  desired  in  engineering  structural  application. 
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PHYSICAL  MODEL:  A  CRACKED  SPECIMEN 

In  order  to  exhibit  the  amount  of  irreversibility  caused  by  physical  damage 
of  the  loaded  material,  the  behavior  of  a  center  cracked  panel,  Figure  1,  under 
unidirectional  tensile  loading  perpendicular  to  the  crack  is  depicted.  Such  a 
specimen  features  a  macro-defect  in  the  form  of  a  symmetrically  located  crack, 
whose  length  is  of  the  same  order  of  magnitude  as  the  specimen  size.  Under  suf¬ 
ficiently  high  loads,  the  macrocrack  would  be  expected  to  increase  in  length. 

For  some  value  of  the  load,  the  process  will  become  unstable  and  result  in  fail¬ 
ure  of  the  panel . 

A  specimen  made  of  a  "ductile"  material  would  be  expected  to  develop  damage 
at  a  lower  scale  level  in  the  form  of  microcracks  and/or  microvoids.  These  de¬ 
fects  may  be  initially  present  in  the  material,  or  may  coalesce  from  the  move¬ 
ment  of  dislocations  and  other  defects  at  an  even  smaller  scale,  i.e.,  the  atomic 
scale.  In  the  presence  of  a  macrocrack,  such  regions  are  known  to  exist  ahead 
of  the  crack  tip  on  both  sides  of  a  two-ended  crack.  The  latter  regions,  common¬ 
ly  identified  as  "plastic  enclaves"  in  the  classical  theory  of  plasticity,  have 
been  the  suDject  of  much  attention  since  the  dimensions  of  such  regions  can  be 
significant  when  compared  to  the  crack  length.  The  microcracks  themselves  will 
produce  microplastic  regions  ahead  of  both  crack  tips,  also.  Their  effect  on 
the  propagation  of  the  macrocrack  is  not  yet  fully  understood. 

In  order  to  develop  a  model  of  the  processes  described  above,  it  is  assumed 
that  the  strain  energy  density,  dW/dV,  per  unit  volume  of  material,  is  the  quan¬ 
tity  that  reflects  material  damage  at  the  micro  level,  as  well  as  the  length  and 
direction  of  macrocrack  propagation  at  the  macro  level.  The  term  "material  dam¬ 
age"  will  refer  to  processes  which  utilize  strain  energy  at  the  continuum  scale 
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Figure  1  -  Center  cracked  panel  specimen 


to  produce  irreversible  changes  in  the  material.  These  changes  will,  in  turn, 
affect  the  lorn1  material  behavior  as  a  result  of  the  formation  of  microcracks, 
microvoido  or  coalescence  of  these  defects.  A  survey  of  the  effective  material 
properties  of  bodies  with  statistically  isotropic  microcrack  arrays  is  discussed 
in  the  first  part  of  this  report  [6]. 

What  is  essential  is  to  establish  a  unique  correspondence  between  each  point 
on  a  load  displacement  curve  while  the  cracked  panel  is  damaged  by  the  build  up 
of  microcracks  in  the  zones  referred  to  as  "macro-plasticity"  in  Figure  1  and 
the  advancement  of  the  macrocrack.  The  results  will  be  informative  for  develop¬ 
ing  constitutive  relations  of  different  crystalline  materials  characterized  by 
the  different  degrees  of  microdamage  and  macrocracking.  In  particular,  it  may 
be  possible  to  distinguish  those  constitutive  coefficients  associated  with  de¬ 
scribing  material  behavior  from  those  characterizing  material  properties.  For 
instance,  t.,e  Young's  modulus  and  Poisson's  ratio  are  the  only  two  material  con¬ 
stants  for  an  isotropic  elastic-plastic  material  while  the  nonlinear  behavior 
is  described  py  the  strain  hardening  exponents  which  are  sensitive  to  change  in 
specimen  sizes. 
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FINITE  ELEMENT  FORMULATION 


All  stress  analysis  used  to  generate  the  load-displacement  curves  were  car¬ 
ried  out  on  center  cracked  panel  specimens  for  materials  with  characteristics 
of  a  low  yield  stress  (a  =  35,000  psi)  structural  steel.  The  load-displacement 
curves  were  calculated  incrementally  in  a  series  of  "stress  analysis-material 
microdamage  assessment-stress  analysi s-macrocrack  extension"  increments  or  cycles. 

The  Axisy...metric/Planar  Elastic  Structures  (APES)  two-dimensional  fracture 

mechanics  and  stress  analysis  finite  element  program  [7]  performed  the  stress 

analyses  portion  of  each  increment.  The  formulation  of  APES  utilizes  quad-12 

elements  which  allow  for  cubic  displacement  fields  and  quadratic  stress  and 

1  /  2 

strain  -Helds  within  each  element.  The  r  displacement  field  in  the  immediate 
vicinity  of  the  crack  tip  is  embedded  in  the  solution  through  the  use  of  1/9  to 
4/9  nodal  spacing  on  the  element  sides  adjacent  to  the  crack  tip.  The  tensile 
loading  of  the  panel  subjects  the  crack  to  a  symmetric  Mode  I  situation,  thus 
requiring  only  one  quarter  of  the  physical  problem.  Figure  2,  to  be  modeled  by 
440  nodes  which  define  77  elements.  The  CDC  6400  computer  used  in  this  study 
typically  required  225  system  seconds  and  132k  of  central  memory  to  complete 
the  stress  analysis  for  a  given  specimen  load  and  damage  distribution. 

The  microdamage  processes  that  occur  in  the  material  manifest  themselves  by 
changes  in  the  material  properties  at  the  continuum  level.  The  finite  element 
formulation  allows  for  this  by  associating  with  each  element  one  of  up  to  25  dif¬ 
ferent  material  property  pairs  (i.e.,  elastic  modulus  E,  and  Poisson's  ratio  v). 
Thus,  a  set  of  properties  (Ep^)  corresponding  to  the  initial  elastic  behavior 
of  the  material,  as  defined  from  an  uniaxial  tensile  test  specimen,  and  up  to  24 
discrete  material  property  pairs  (E^v^  for  i  =  2,3,4, ...  ,24,  different  instances 
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in  the  damaged  material's  developing  history. 

The  extension  of  the  macrocrack  in  Figure  1  is  analyzed  by  convecting  the 
grid  boundaries  in  the  finite  element  analysis  for  each  increment  of  macrocrack 
growth.  This  is  done  so  that  the  crack  tip  and  surrounding  region  are  not  dis¬ 
torted  when  the  crack  grows.  Otherwise,  distortion  of  the  finite  elements  would 
interfere  with  the  grid  pattern  caused  by  material  damage.  Convection  of  the 
boundaries  confines  any  element  distortion  to  the  boundary  elements  which  are 
less  likely  to  sustain  microdamage  at  the  loads  considered. 
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MICRODAMAGE  MODEL 


The  damage  model  used  for  the  center  cracked  panel  specimen  was  based  on  24 
discrete  damage  levels  of  the  virgin  material.  (Recall  that  25  different  mate¬ 
rials  can  be  incorporated  into  the  finite  element  formulation).  This  model  will 
be  referred  to  as  the  Multiple  Damage  Level  Model.  The  24  damage  levels  were 
defined  on  the  basis  of  equal  multiple  elastic  secant  modulus  reductions.  Figure 
3.  The  24th  is  defined  from  the  ultimate  stress  and  strain,  au  and  eu,  of  the 
true  stress- true  strain  curve.  Since  each  of  the  24  damage  levels  is  attained 
at  a  different  value  of  strain  energy  density,  the  true  stress-true  strain  curve 
must  be  explicitly  given.  Figure  4.  For  this  analysis,  the  form  used  was 


*  =  E] e  *  <  °y 


a  =  Ae  a  >_ 


For  both  materials,  the  stress-strain  relations  are 


(1) 


j  -  30 .0x1 0^  c  a  <  cfy 
a  =  190500  c0,251  a  >  cy 


(2) 


which  is  based  on  a  structural  steel.  For  the  onset  of  damage  in  the  material, 
the  strain  energy  density  in  the  ith  element  must  be  greater  than  the  strain  en¬ 
ergy  density  at  the  yield  stress 


(—) 

lAVj. 


/  AWs 
(AV'. 


(3) 
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where  the  strain  energy  density  for  the  case  of  plane  strain  is  given  by 


AW 

AV  "  2E 


(4) 


For  the  k^n  level  of  damage  in  the  material,  characterized  by  an  elastic  modulus 
of  Ek,  the  strain  energy  density  in  the  ith  element  must  be 


/ AW ,  / AW  \  /AW< 

<  (A\T>  1  AV" 


k+l 


(5) 


where,  for  the  form  of  the  given  stress-strain  relation 


(w)  =  /  -de  =  /  adz  +  j  ode 


lAV 


k  o 


/ AW\  ,  A  /  n+1  n-1 , 

=  n+T  ^£k  -  ey  ) 


1  Jy  ,  A  ,  n+1  n-1 , 

=  2  4  +  HTT(ek  -  ey  } 


(6) 


The  load  P  is  the  independent  variable  for  this  model,  since  all  elements 
have  the  potential  to  sustain  damage  at  a  given  load.  The  increment  of  the  load 
P  can  be  arbitrarily  chosen  within  the  limits  of  upper  and  lower  bounds  as  de¬ 
scribed  in  [6]. 

The  concept  of  the  mean  damage  level  d  and  damage  center  ,y^)  have  been 
introduced  in  order  to  quantify  the  damage  zone  in  terms  of  a  scalar  magnitude 
and  its  position  in  the  panel  specimen.  The  mean  damage  level  d  is  defined  on 
the  basis  of  the  discrete  reductions  in  elastic  modulus.  Defining  the  damage  in 
the  ith  element  as  the  fractional  reduction  in  the  elastic  modulus  by 
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(7) 


d,  = 


ErEk 


1  E, 


The  mean  damage  level  d  is  calculated  as 


d  = 


T7  ?divi 


(8) 


Elements  which  have  not  reached  the  (dW/dV)  damage  threshold,  of  course,  have 
values  of  d^  equal  to  zero.  The  coordinates  (x^y^)  the  damage  center  are  de¬ 
fined  on  the  basis  of  the  first  moment  of  the  damage  distribution 


x,  = 


x -d ■ V  • 


d  Sd • V •  7  i  i  i 
l  i  1  1 


yd  =  Id  .V  .  |yidiVi  1  =  1 ,2’*  “  ,nc 

l  1  1 


(9) 


where  nd  is  the  number  of  elements  having  damage  at  the  applied  load.  This  is 
simply  the  discretized  form  of  first  moment  integrals  with  the  density  function 
taken  to  be  the  damage  level.  The  and  y.  coordinates  are  the  coordinates  of 
the  ith  element  centroid. 


The  locus  of  the  damage  center  coordinates  at  each  load  serves  as  an  indica¬ 
tion  of  how  the  material  damage,  as  defined  by  the  models  incorporated  here,  is 
being  distributed  in  the  center  cracked  panel.  The  crack  tip  is  the  origin  of 
the  locus  of  points,  and  the  subsequent  curve  is  indicative  as  to  whether  the 
lower  levels  of  damage  at  the  periphery  of  the  damage  zone,  or  the  higher  levels 
nearer  to  the  crack  tip  are  predominating  the  damage  accumulation  process  at  a 
given  load  increment.  Also,  the  effect  of  the  free  boundary  ahead  of  the  crack 
can  be  evaluated  as  to  its  ability  to  influence  the  damage  zone's  direction  of 
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propagation.  Typically,  one  would  expect  the  damage  zone  boundary  to  propagate 
more  toward  the  free  surface  as  the  zone  si2e  becomes  significant  with  respect 
to  the  uncracked  ligament  length. 
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STRAIN  ENERGY  DENSITY  THEORY  OF  MACROCRACK  PROPAGATION 


Strain  Energy  Density  Theory,  as  proposed  by  Si h  [4,5],  makes  use  of  the  en¬ 
ergy  state  and  material  properties  ahead  of  the  crack  for  predicting  the  direc¬ 
tion  and  increment  of  crack  growth.  Physically,  the  material  behavior  in  the 
immediate  vicinity  of  the  crack  tip  is  too  complex  to  be  described  by  a  contin¬ 
uum  mechanics  analysis  because  of  the  influence  of  material  microstructure. 

This  necessitates  the  investigation  to  be  centered  on  elements  at  a  finite  dis¬ 
tance  away  from  the  crack  tip  such  that  the  material  properties  can  be  described 
by  parameters  at  the  continuum  or  macroscopic  scale  level. 

The  strain  energy  density  field  in  the  neighborhood  of  the  crack  tip  has  the 

form 


[S(a,e;material  properties)]  +  nonsingular  terms  (10) 

where  r  is  the  radial  distance  from  the  crack  tip,  a  is  the  half  crack  length, 
e  is  the  angular  position  coordinate,  and  the  function  S  is  the  strain  energy 
density  factor  which  is  the  coefficient  of  the  r"^  singular  energy  density  field. 
In  the  case  of  linear,  homogeneous  and  isotropic  elasticity,  the  strain  energy 
density  factor  can  be  expressed  in  terms  of  the  stress  intensity  factors  k-j ,  k2 
and  k^  as 

S  =  ank|  +  a12k]k2  +  a22k^  +  a-^k^  (11) 

where  the  a—  (i,j  =  1,2,3)  are  functions  of  the  elastic  shear  modulus,  Poisson's 
ratio,  and  o.  The  k^-factors  take  the  form 
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k  .  =  c.  av'a 

J  J 


(12) 


where  the  c.  are  constants  that  depend  on  the  solid  geometry  and  a  is  the  load- 

J 

ing.  In  the  case  of  yielding,  no  realistic  analytical  crack  tip  stress  or 
★ 

strain  expressions  are  available.  Because  of  the  nonhomogeneous  nature  of  the 
elastic-plastic  material  behavior  near  the  crack  tip,  a  single  parameter  repre¬ 
sentation  such  as  a  stress  intensity  factor  does  not  exist.  Numerical  calcula¬ 
tion  is  required  for  analyzing  the  state  of  affairs  near  the  crack  tip. 

A  quasi -linear  approach  is  taken  such  that  each  segment  of  crack  growth  ac¬ 
companied  by  microdamage  is  analyzed  using  the  l//r  singularity  solution.  The 
global  response  consisting  of  all  the  increments  of  crack  growth  is  nonlinear 
because  the  material  will  not  be  damaged  uniformly  around  the  crack.  This  ef¬ 
fect  is  assessed  quantitatively  by  a  numerical  procedure  calculating  the  local 
strain  energy  density  field. 

The  fundamental  hypotheses  of  the  strain  energy  density  theory  applied  to 
both  fracture  and  yielding  may  be  stated  as  follows: 

(I)  Yielding  is  assumed  to  coincide  with  the  location  of  relative  maximum 
of  AW/AV,  i.e.,  ( AW/ aV )  and  to  occur  when  (aW/aV)  reaches  a 

inaX  iTlaX 

critical  value  (aW/aV)  . 


if 

Available  solutions  [8]  assume  that  the  material  yields  uniformly  around  the 
crack  tip,  a  situation  that  deviates  far  from  reality. 
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(II)  Fracture  is  assumed  to  be  determined  by  the  location  of  relative 


rr  ..imum  of  (aW/aV),  i.e.,  (AW/AV)m.jn  and  to  occur  when  (AW/AV)m^n 
reaches  a  critical  value*  (aW/aV)c. 


★ 

As  mentioned  earlier,  if  the  material  along  the  prospective  fracture  path  is 
also  yielded  or^damaged  microscopically,  then  the  critical  value  must  be  modi¬ 
fied  as  (aW/aV)c  depending  on  the  severity  of  material  damage. 
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STRESS  AND  FAILURE  ANALYSIS 


The  development  of  damage  zones  and  the  propagation  of  a  macrocrack  are 
two  irreversible  damage  processes  which  must  be  included  to  model  the  behavior 
of  a  real  ductile  material.  Since  material  damage  is  a  load  history  dependent 
process,  stress  and  failure  analysis  must  be  performed  in  tandem  for  each  incre¬ 
ment  of  loading.  Referring  to  the  cracked  panel  in  Figure  1,  the  load  is  in¬ 
creased  until  the  strain  energy  density  function  (aW/aV)  for  elements  in  the  vi¬ 
cinity  of  the  crack  reaches  a  level  corresponding  to  uniaxial  tensile  yield,  i.e., 
(aW/aV)^.  Thus,  a  contour  of  (aW/aV)^  =  const,  will  be  developed  and  enclose  a 
larger  region  for  each  increment  of  loading.  Macrocrack  propagation  is  assumed 
to  start  when  (aW/aV)  reaches  a  critical  value  (aW/aV)  .  As  an  example,  if  the 
macrocrack  happens  to  advance  into  the  undamaged  material,  the  onset  of  macrocrack 
growth  is 


* 


On  the  other  hand,  had  damage  already  developed  along  the  prospective  fracture 
path,  say  by  the  amount  (AW/AV)^,  then,  the  condition  of  macrocrack  growth  be¬ 
comes 


★ 

/ AW,  _  , AW,  ,AW, 

'ait  =  ' av  '  'Air  . 

c  c  d 


(14) 


Equation  (14)  can  thus  be  repeated  for  each  successive  increment  of  crack  growth. 

The  relation  of  (aW/aV)^  to  (AW/AV)  for  the  tensile  specimen  is  illustrated 
schematically  in  Figure  5.  The  area  OYA  represents  the  energy  (aW/aV)^  which 
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Figure  5  -  Strain  energy  density  definitions  from  uniaxial 
true  stress-true  strain  curve 

has  been  dissipated  by  the  microdamaqe  process  due  to  loading  along  the  stress 

* 

strain  curve  up  to  point  A.  The  relative  critical  strain  energy  density  (aW/aV)c 
is  the  required  amount  of  energy  necessary  for  the  critical  strain  energy  density 
(tW/cV)c  level  to  be  reached  in  the  material. 

If  the  specimen  were  unloaded,  the  pseudo-linear  elastic  nature  of  the  as¬ 
sumed  material  model  would  cause  unloading  to  proceed  along  the  line  AO  back  to 
the  origin.  The  area  OAB  is  then  the  local  recoverable  elastic  energy  ( aW/^V ) e . 

If  loading  were  to  continue  further  along  curve  AU,  the  additional  energy  required 
to  reach  (iW/AV)c  is  (aW/aV)  in  view  of  equation  (14).  Thus,  the  relative  criti¬ 
cal  strain  energy  density  can  also  be  expressed  as  the  sum  of  the  local  recover¬ 
able  and  additional  energy  densities,  i.e.. 
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(15) 


(AW  MV)*  =  (AW /AV)  +  (AW/AV)* 


To  reiterate,  microdamage  is  assumed  to  occur  in  the  region  local  to  the 

macrocrack  front  when  a  material  element  modeled  by  finite  elements  exceeds  a 

certain  limit  (aW/aV)  .  Based  on  the  strain  energy  density  criterion,  macrocrack 

★  •  .  _ 

growth  begins  when  the  relative  minimum  of  the  strain  energy  density  function 
reaches  the  condition  in  equation  (14).  The  stable  crack  growth  process  con¬ 
tinues  until  the  condition  of  global  instability 


★ 


★ 


(16) 


★ 

is  reached.  The  parameter  rc  represents  the  last  crack  growth  increment  that 

★ 

corresponds  to  global  instability  and  Sc  may  be  interpreted  as  the  fracture 
toughness  of  the  material. 


The  load-displacement  curves  of  the  cracked  panel  are  produced  by  successive 
applications  of  the  APES  finite  element  program  for  each  load  increment.  This 
procedure  incorporates  the  stress  and  failure  analysis  involving  microdamage 
and  crack  propagation  either  separately  or  in  combination  with  the  method  de¬ 
scribed  above.  A  flow  chart  describing  the  step-by-step  procedure  is  shown  sche¬ 
matically  in  Figure  6. 


Even  when  the  material  surrounding  the  crack  tip  is  completely  yielded,  the  di¬ 
rection  of  crack  initiation  is  still  determined  by  (aW/aV)  •  calculated  in  the 
yielded  material . 
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LOAD-DISPLACEMENT  RESPONSE 


The  combined  effects  of  microdamage  and  crack  propagation  on  the  incremental 
load-displacement  curves  are  shown  in  Figure  7  for  two  materials  with  different 
(aW/aV)c  values.  The  nonlinear  response  owing  to  damage  is  seen  in  terms  of  the 
departure  of  the  curves  from  the  linear  response.  According  to  the  quasi-1 inear 
model,  the  difference  of  the  nonlinear  and  linear  responses  of  the  panel  with  a 
material  (aW/aV)  value  of  466.7  lb-in/in3  are  seen  to  be  approximately  60"  of 
those  for  the  panel  with  a  material  (aW/aV)c  value  of  64.2  lb-in/in3. 

In  Figures  8  and  9,  the  separated  and  combined  effects  of  microdamage  and 
macrocrack  propagation  are  shown  for  both  materials.  The  nonlinear  response  of 
both  panels  due  to  the  separate  and  combined  effects  of  material  damage  at  both 
scale  leva's  are  shown.  As  would  be  expected,  the  individual  contributions  of 
microdamage  and  macrocrack  propagation  are  less  than  the  combined  effect.  Of 
particular  interest  are  the  proportions  of  the  individual  effects  for  the  two  ma 
terial s . 

In  the  case  of  the  less  tough  material,  with  (aw'/aV)c  =  64.2  lb-in/in3,  the 
global  displacements  are  greater  for  macrocrack  propagation  alone  with  no  micro¬ 
damage  as  compared  to  the  case  of  microdamage  without  crack  propagation.  Thus, 
microdamage  appears  to  have  less  influence  on  the  global  panel  stiffness  than 
that  due  to  crack  propagation.  For  the  high  toughness  material  with  (aW/aV)c 
=  466.7  lb-in/in3,  the  reverse  trend  is  observed.  The  microdamage  effects,  when 
taken  alone,  are  seen  to  produce  a  greater  reduction  in  global  panel  stiffness 
than  macrocrack  propagation. 

The  foregoing  load-displacement  curves  were  generated  from  (aW/aV)c  values 
for  the  undamaged  portion  of  the  material  as  defined  from  the  area  under  the  uni 
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axial  true  stress-true  strain  curves.  While  the  microdamage  process  is  a  phys¬ 
ically  irreversible  process,  crack  propagation  is  a  sequence  of  material  separa¬ 
tion  events  contingent  on  the  current  damage  distribution  and  relative  critical 

★ 

strain  energy  density  function  (aW/aV)c  throughout  the  material.  Clearly,  in  a 
material  with  comparatively  lower  (aW/aV)c,  the  damage  process  takes  place  over 
a  narrower  range  of  strain  energy  density  values.  In  a  material  with  higher 
(aW/aV)  ,  this  range  of  energies  is  greater,  allowing  for  higher  levels  of  dam¬ 
age  prior  to  material  degradation  coincident  with  fracture. 

The  variations  of  global  secant  stiffnesses  with  load  for  microdamage  and 
macrocrack  propagation  are  shown  in  Figures  10  and  11.  The  comparative  reduction 
in  stiffness  for  the  individual  and  combined  effects  of  microdamage  and  crack 
propagation  are  again  apparent.  For  a  (aW/aV)c  value  of  64.2  lb-in/in3,  and  at 
the  maximum  load  considered,  the  effect  of  microdamaqe  is  seen  to  reduce  the  pan¬ 
el  stiffness  by  -10"  while  for  the  case  of  crack  propagation,  the  stiffness  is 
reduced  by  -20".  The  combined  effects  are  greater  than  either  taken  individually, 
reducing  the  stiffness  by  ~33/o.  For  a  (AW/AV)C  value  of  466.7  lb-in/in3,  the 
stiffness  reductions  for  microdamage,  propagation  and  combined  effects  are  ap¬ 
proximately  16,  6  and  23%,  respecti vely,  at  the  maximum  load  considered.  The 
relative  influence  of  microdamage  is  seen  to  be  greater  than  the  effect  of  crack 
propagation,  as  discussed  above  on  the  basis  of  the  load-displacement  curves. 

In  Figures  12  and  13,  the  propagation  of  the  macrocrack  is  shown  in  the  ab¬ 
sence  of  any  microdamage.  (For  this  reason,  only  a  portion  of  the  finite  ele¬ 
ment  grid  is  used  to  illustrate  the  propagation  of  the  crack).  This  case  corre¬ 
sponds  to  an  ideally  brittle  material  which  fractures  along  flat  surfaces.  As 
would  be  expected,  the  crack  propagates  further  for  the  material  with  a  lower 
(aW/aV)c  value.  Stable  crack  growth  has  been  assumed  in  this  analysis,  and  in 
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all  subsequent  discussions  involving  crack  growth. 

The  microdamage  zones  in  the  cracked  specimens,  as  characterized  by  the 
original  material  elastic  modulus  of  30.0  x  10^  psi  and  twenty-four  damaged  ma¬ 
terial  states  for  the  finite  element  discretization,  are  shown  in  Figures  14  and 
15.  Recall  that  material  no.  25  is  the  highest  level  of  damage  in  the  discreti¬ 
zation,  with  a  tensile  specimen  secant  modulus  of  18.0  x  10^  psi  for  a  (aW/aV)c 
value  of  64.2  lb-in/in3,  and  6.0  x  10^  psi  for  a  (aW/aV)c  of  466.7  lb-in/in3. 

All  elements  without  a  number  are  assumed  to  have  not  reached  the  yield  condition, 
based  on  the  average  strain  energy  density  taken  over  the  element,  and  hence  are 
material  no.  1. 

In  Figures  16  and  17,  the  extent  of  microdamage  is  shown  when  crack  propaga¬ 
tion  is  also  taken  into  consideration.  No  major  differences  in  the  damage  zones 
for  each  material  can  be  seen  between  a  propagating  and  non-propagating  macrocrack. 
There  is  some  enhancement  in  the  damage  levels  of  the  combined  damage  and  crack 
propagation  cases,  as  would  be  expected  due  to  increasing  macrocrack  length.  The 
damage  zones  in  both  cases,  for  both  materials,  are  seen  to  be  similar  in  the 
shape  of  "plastic  enclaves"  predicted  by  the  von  Mises  criterion  (which,  as  dis¬ 
cussed  above,  is  not  based  on  microdamage,  but  rather  traditional  continuum  me¬ 
chanics),  and  observed  in  fracture  specimens. 

The  crack  growth  increments  shown  above  were  determined  according  to  the 
strain  energy  theory  as  discussed  previously.  The  strain  energy  density  fields 
ahead  of  the  crack  tip,  with  and  without  the  effects  of  microdamage,  at  each 
load  increment  are  shown  in  Figures  18  and  19.  The  predominant  characteri sties 
of  these  families  of  curves  are  the  increasing  intensity  of  the  energy  field  for 
increasing  crack  lengths  and  the  decreasing  value  of  the  relative  critical  strain 
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energy  (aW/aV)c  as  the  damage  level  ahead  of  the  crack  tip  increases.  Due  to 
the  discretization  of  the  panel  domain  required  by  the  finite  element  method, 
the  (aW/aV)c  value  is  treated  as  a  locally  constant  quantity  in  each  individual 
element.  This  is,  of  course,  similar  to  the  way  the  elastic  moduli  of  the  in¬ 
dividual  elements  are  treated. 

The  average  secant  modulus  in  the  microdamage  zone  is  plotted  as  a  function 
of  load  in  Figure  20.  For  both  materials,  with  and  without  macrocrack  propaga- 


Figure  20  -  Average  secant  modulus  in  microdamage  zone 
versus  applied  load 
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tion,  the  average  secant  modulus  within  the  damage  zone  is  seen  to  be  between 
83  and  88  percent  of  the  undamaged  value.  This  is  interesting  in  light  of  the 
lowest  values  of  the  secant  moduli  characteristic  of  the  two  materials  considered. 
These  values  are  0.6  and  0.2  times  the  undamaged  modulus.  The  implication  is 
that  the  effect  of  the  higher  levels  of  damage  (particularly  near  the  crack  tip) 
on  the  average  modulus  are  balanced  by  the  lower  levels  of  damage  near  the  ex¬ 
panding  periphery  of  the  damage  zone  itself. 

The  damage  center  loci  (Figure  21)  for  each  material  also  show  similar  be- 

5 

havior.  For  loads  less  than  1.8  x  10  pounds,  the  loci  follows  the  direction 
0Q  =  cos”^  (l-2v)  of  maximum  strain  energy  density  in  the  neighborhood  of  the 
crack  tip,  as  discussed  by  Sih  [4]  for  a  linear  elastic  body.  The  departure  of 

5 

the  loci  from  this  direction  occurs  at  applied  loads  greater  than  1.8  x  10 
pounds  in  the  direction  of  the  free  surface.  This  is  seen  to  be  independent  of 
the  value  of  (aW/AV)c,  or  whether  macrocrack  propagation  is  included  as  a  dissi¬ 
pative  mechanism.  For  the  panel  with  a  higher  value  of  (aW/aV)c,  the  movement 
of  the  loci  away  from  the  9q  direction  is  seen  to  be  greater  than  for  the  panel 
with  the  lower  value  of  (aW/aV)c.  Thus,  the  increased  dissipative  capability 
characteristic  of  the  higher  (aW/AV)c  material  extends  the  damage  zone  toward  the 
free  surface.  This  relative  movement  of  the  damage  zone  may  be  related  to  the 
onset  of  stable  and  unstable  macrocrack  growth.  It  may  be  produced  as  the  re¬ 
sult  of  the  damage  zone  size  and  geometry  itself,  the  influence  of  the  free  sur¬ 
face  ahead  of  the  crack  tip,  or  a  combination  of  the  above. 
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CONCLUDING  REMARKS 


The  global  behavior  of  a  center  cracked  specimen  has  been  modeled  by  the 
finite  element  method  while  allowing  for  material  damage  at  each  load  increment. 
Both  microdamage  (in  the  form  of  isotropic  distributions  of  microcracks)  and 
macrocrack  propagation  have  been  modeled  by  a  consistent  application  of  the 
strain  energy  density  criterion.  The  thresholds  of  material  damage  are  referred 
to  data  from  the  true  stress-true  strain  curves  of  uniaxial  tension  tests. 

The  nonlinear  behavior  of  the  load-displacement  curves  for  both  separate  and 
combined  application  of  the  microdamage  and  macrocrack  propagation  effects  have 
been  incrementally  generated  for  two  different  materials  as  characterized  by 
their  critical  strain  energy  density.  Both  the  load-displacement  curves  and 
secant  stiffness  curves  are  material,  specimen  geometry  and  load  configuration 
dependent.  In  this  analysis,  only  the  material  properties  of  the  cracked  speci¬ 
mens  have  been  changed  in  accordance  with  the  severity  of  the  damage.  Conceptu¬ 
ally,  this  approach  could  be  applied  to  develop  constitutive  relations  for  mate¬ 
rials  in  which  damage  due  to  micro  and  macrocracking  can  be  treated.  In  particu¬ 
lar,  the  procedure  for  a  quantitative  assessment  of  material  damage  in  terms  of 
load-displacement  response  has  been  established. 

The  results  of  the  analyses  show  how  macrocrack  propagation  can  dominate  the 
global  nonlinearity  when  the  critical  strain  energy  density  level  is  relatively 
low.  Conversely,  for  the  material  with  a  higher  critical  strain  energy  density, 
the  microdamage  effect  becomes  more  influential.  These  two  observations  are  con¬ 
sistent  with  the  generally  observed  results  for  materials  that  behave  in  a  pri¬ 
marily  brittle  and  primarily  ductile  fashion. 
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On  the  basis  of  the  information  obtained  from  the  present  analysis,  coupled 
with  additional  and  more  refined  damage  models,  appropriate  tensile  tests  may  be 
designed  to  reflect  the  damage  processes  which  are  actually  occurring  in  the 
structural  component.  The  appropriate  choice  of  variables,  which  will  be  capa¬ 
ble  of  describing  the  s^ate  of  material  damage,  can  then  be  incorporated  into 
viable  constitutive  relations.  These  relations  are  useful  for  predicting  the 
global  behavior  of  structural  members  from  material  damage  at  the  continuum  scale. 
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